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Abstract 

We show a function that fits well the probability density of return times between two consecutive 
visits of a chaotic trajectory to finite size regions in phase space. It deviates from the exponential 
statistics by a small power-law term, a term that represents the deterministic manifestation of the 
dynamics. We also show how one can quickly and easily estimate the Kolmogorov-Sinai entropy 
and the short-term correlation function by realizing observations of high probable returns. Our 
analyzes are performed numerically in the Henon map and experimentally in a Chua's circuit. 
Finally, we discuss how our approach can be used to treat data coming from experimental complex 
systems and for technological applications. 
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Observing how long a dynamical system takes to return to some state is one 
of the simplest ways to model and quantify its dynamics from data series. In this 
work, we describe a simple way to extract some relevant invariant quantities of 
a chaotic system by using recurrence times, in particular Poincare recurrences 
that measure the time interval for a system to return to a configuration close to 
its initial state. Part of this work is dedicated to apply the theoretical results 
proposed by [Pinto et al. arXiv:0908.4575] to calculate the Kolmogorov-Sinai 
entropy and the decay of correlation of the experimental Chua's circuit when 
the returns are measured in "large" size regions in phase space. Another part is 
dedicated to study how the deviation (from the exponential form) of the density 
of the first Poincare returns can be used to detect deterministic manifestations in 
chaotic systems. Finally, we discuss how our approach can be used to treat data 
coming from experimental complex systems and for technological applications. 



I. INTRODUCTION 

Chaotic systems have simultaneously a stochastic [1] and a deterministic dynamical char- 
acters {2]. Single trajectories are predictable (deterministic) for a short-term evolution and 
unpredictable for a long-term evolution (stochastic). While the stochastic character is asso- 
ciated with an exponential decay of correlations and information about the actual state is 
rapidly lost as the trajectory evolves, the deterministic character is associated to a power-law 
decay of correlations and information about the actual state is lost as the system is evolved, 
but in a rather slower fashion. 

A relevant question that arises when dealing with data coming from complex systems is 
whether manifestations of chaotic behavior can be detected in the data [3] so that one can 
construct deterministic models. One wishes to come up with a dynamical description of the 
data, but due to the sensitivity dependence on initial conditions of chaotic systems one is 
prompt to adopt a probabilistic approach in order to reveal the underlying dynamics of the 
system from statistical averages 4j. A promising tool of analysis is provided by the statistics 
of the Poincare recurrence time (PRT) 5] which measures the time interval between two 
consecutive visits of a trajectory to a finite size region in phase space. 

Many relevant quantifiers of low-dimensional chaotic systems can be obtained by the 



2 



statistical properties of the PRTs. The purpose of the present work is to apply some results 
from Ref. 6| and to propose other theoretical approaches to easily identify deterministic 
and stochastic manifestations in dissipative strongly chaotic systems by using the PRTs of 
chaotic trajectories to regions of finite size (7|, considering only short return times. These 
later two conditions constraining our analyzes are devoted to suitably apply our approaches 
to realistic physical situations: the resolution to measure returns as well as the time frame 
to realize the experiment is finite. 

By chaotic systems, we mean non-uniformly hyperbolic systems. We focus the analyzes 
on the Henon map and on an experiment, the Chua's circuit, but we also present results 
coming from the logistic map. 

We first present the theoretical framework to be used in Sees. ILTjIHfIV[ and |V] Then, 
in Sec. IVIt we apply this approach to analyze the Henon map and the experimental Chua's 
circuit . We also work with the Logistic map in Appendix |A] 

In Sec. HIl we present a continuous function, pp that fits well the probability distribution 
for the PRTs to regions of finite size. The many parameters of pp are theoretically estimated 
in Sec. IIHI Further, we show how to use these parameters to quantify how stochastic or 
deterministic the considered system is, under the perspective of the PRTs. Concerning 
the coefficient of the power-law term (responsible for the deviation of pp from the Poisson 
distribution), a first approximation furnishes that it is inversely proportional to the average 
return time. The longer (shorter) is the average return time, the smaller (larger) the power- 
law term, and the more stochastic (more deterministic) the PRTs. 

In Sec. IIVI we succinctly describe how to calculate the Kolmogorov-Sinai (KS) entropy, 
denoted by H KS , in terms of the number of unstable periodic orbits (UPOs) and of the 
frequency with which the PRT happen. Then, we explain how Hks can be calculated in 
typical physical situations, when the information about the UPOs is unknown and the only 
available information is the frequency with which a high probable PRT with short length 
happens. Notice that this frequency can be easily measured even if only a few return times 
are observed because faster returns are more probable. 

Alternative methods to calculate the correlation entropy, K 2 , a lower bound of Hxs-, 
and to calculate H KS from time series were proposed in Refs. [8|, |9j. In Ref. [8| K 2 is 



calculated from the correlation decay and in Ref. 



9| by the determination of a generating 



partition of phase space that preserves the value of the entropy. But while the method in 
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Ref. [8j unavoidably suffers from the same difficulties found in the proper calculation of 
the fractal dimensions from data sets, the method in Ref. ji| requires the knowledge of the 
generating partitions, information that is not trivial to be extracted from complex data. 
The advantage of the theoretical approach used is its simplicity. As one can see in Eq. (|19[) . 
the only information required is the frequency with which high probable PRTs happen in 
some regions of the phase space. 

In Sec. |V] we show how the PRTs can be used to calculate the correlation function for 
short-term returns and an upper bound for long-term returns, a function that indicates how 
much the future returns are related to the past returns to a finite size region. 

Other methods to calculate the correlation function from the PRTs were proposed in 



Refs. 



10 



One of the advantages of the proposed theoretical approach is that we could 



show, in a very trivial way, that if the distribution of the PRTs has an exponential decay 
;hen the correlation function also decays exponentially, a result rigorously demonstrated in 

HI- 

Finally, in Sec. IVII1 we present our conclusions and discuss how this approach can be 
used in an integrated way to characterize experimental data coming from complex systems. 



II. THE PROBABILITY DISTRIBUTION OF THE PRTS 



In recent works 
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13], it was shown that the fitted probability density function Pf(t, B), 
of a series of Poincare return times r n {n — 1, . . . , L) to a box B of equal finite sides 2e 



14] . of typical trajectories in a non-uniform hyperbolic attractor on TZ 2 , deviates from the 



exponential law pp(r, B) = p(B)e~ TfJj( - B ^ if the box B is placed on some special region of the 
phase space, where p{B) is the probability measure inside it. This quantity measures the 
frequency with which a chaotic trajectory visits B. 

The function ppij, B) describes well the discrete probability distribution function (PDF), 
denoted by p(r, B), in Axiom-A systems and some classes of non- uniformly hyperbolic sys- 



tems with exponential decay of correlations 
arbitrarily small intervals (e — > 0). In Ref. 
have p(r, B) — > p' F (r, B) with 



=chaos with strong mixing properties) 15] , for 
13] , it was hypothesized that for r > r^ p n , we 



p' F (r,B) = (3e 



(-or) 



(1) 



where r > T™ p n represents the PRT for which the distribution p(r, B) becomes approxi- 
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mately continuous and /3 = (3r + a. It was assumed that f3^e^~ aT ^ describes an effective 
PDF associated mostly with the return of trajectories along the unstable manifold of UPOs 
outside the box B. It is the result of the non-hyperbolic nature of the dynamics inside the 
box, whose probability measure suffers the influence of non-local UPOs. The term ae'" aT ' 
describes the hyperbolic nature of the dynamics inside the box, the return of trajectories to 
the box associated with the local dynamics provided by the UPOs inside the interval. 

The coefficients of p' F {r,B) are obtained by the least square fitting method which mini- 
mizes the error between p' F (r, B) and p(r, B): 

E(p' F - p) = (p' F - pf . (2) 



In addition, assuming that Eq. flTJ describes well p(r,B), then the following equations 
should be simultaneously satisfied 

—min 

p{r,B)+ / P ' F {r,B)dT = 1 (3) 



nun 
UPO 



■nun 
T UPO 



p(r,B) x r + / Mr,B)Tdr = (r) (4) 

7~m in UPO 

by considering that p(r, B) can be broken in two terms, one that describes its discrete 
nature, the probability of finding a PRT of length r < T^ P n 0) and another that describes its 
continuous nature. 

Denoting r m j n as to be the PRT with the minimum length [min (r n )] and (r) = 
l / L Y.n=i T m and assuming that r min = r^p then a = ((r) - r min ) _1 and (3 R = 
a [ e (aT rara ) _ ]^ f Qr sufficiently small e so that the terms that contain max(r n ), regarded 
as r max , can be neglected. So, for finite (r), (3 > a. However, Eq. flTJ does not seem to com- 
pletely capture the nature of p(r, B) for finite size intervals. Satisfying conditions flHJ) and 
fllj) do not necessarily minimizes the error in fl2J) , and the contrary also does not apply. Such 
disagreement becomes stronger in boxes centered close to homoclinic tangencies. By fitting 
p(r, B) by a function of the type in Eq. flTJ) might produce (3 < a, which disagrees with 
the inequality (3 > a that validates Eq. flTJ. For the here considered dynamical systems for 
parameters far away from intermittent behavior, and data coming from plasma turbulence 
and stock market [3], one finds often that a = (3, a consequence of p(r,B) be greater than 
p' F (r,B), for t < (r), and that p(r,B) < p' F (r,B), for r > (r). These facts suggest that 



for systems with a dynamics similar to the Henon map, the small term (3r considered to 



be constant in Ref. 



13( is in fact a function that we hypothesize to be a power-law with 



respect to t' . That lead us to 

Pf (t',B)=((3* + c 

where 



t' 



(Te) 



e ( " aT,) (5) 



'UPO 



1 (6) 



(r e ) = / p(r\B)r'dr' (7) 



mm 
UPO 



and 7 is a positive small value. Notice that the use of r' and (r e ) is only an artifact to 
simplify our further approximations and also to simplify any possible nonlinear fitting that 
we might make concerning the real PDF. 

The power-law term /3r = c[l — (-p-y) 7 ] shows that the considered systems for which this 
distribution holds have a return time distribution whose decay is not characterized by either a 
power-law or exponential decay but by both. It appears as a combined effect of the finiteness 
of the box, the expanding factor nearby low-periodic UPOs along the unstable direction, and 
the existence of homoclinic tangencies, which gives the non-uniformly hyperbolic character 
of the Henon map and the Chua's circuit. 

A large value of 7 (which implies in a large /3r) indicates that there is a large contribution 
to the measure inside the box due to UPOs that are outside the box. Two things contribute 
for a large value of 7: the size and location of the box. The larger the size and the level of 
non-hyperbolicity of the box are, the larger the value of 7 is. We find larger values of 7 when 
we measure PRTs in boxes placed close to homoclinic tangencies. At such a case, trajectory 
points that have returned once to the box after P iterations keep consecutively returning 
to the box after P iterations. The trajectory is no longer under the influence of the linear 
expanding character of a period-P UPO (along the direction of the unstable eigenvector) 
but under the influence of the non-linear character of the unstable manifold of the UPO 
located outside the box. 

A sketch of these ideas is depicted in Fig. [TJ In (A), we represent a box that is centered 
on an UPO (in the figure represented by a full circle and the letter O). The box is sufficiently 
large so that the escape of trajectories from the box which are associated with this UPO 
happens no longer along the unstable eigenvector (E u ) but along the unstable manifold 



(A) 



(B) 



W u 



9, 


/ 


/ 


\ 




FIG. 1: A sketch of the structure of the manifolds and eigenvectors of UPOs (represented by full 
circles and the letter O) inside (A) or outside (B) a finite size box, where the PRTs are being 
measured. In (A), the box is far away from homoclinic tangencies, and in (B) the box is in the 
neighborhood of a homoclinic tangency (represented by the full square and the letter T). 



(W u ). The smaller is the period of an UPO, the smaller is the eigenvalue of the unstable 
eigenvector. So, the escape of the trajectory along the unstable manifold is usually expected 
to happen for a low-period UPO. In (B), we represent a box centered in a homoclinic 
tangent (in the figure represented by a full square and the letter T). In the neighborhood of 
homoclinic tangencies exist low-period UPOs whose unstable and stable manifold happens 
to be almost parallel, forming a trapping region. In these trapping regions, the trajectory 
returns to the box along W u and eventually escapes the box, but it may return to it along 
the stable manifold W s of the UPO outside the box. For short returns, these confined 
trajectories contribute positively to p(r, £>), increasing p{B). 

We can classify regions B by the way returns happen inside it: type I are the regions 
whose returns are associated with higher-period UPOs, here conveniently regarded as the 
hyperbolic ones, type II and III are the ones associated with lower-period UPOs, here 
conveniently regarded as the non-hyperbolic regions. In these regions, the deviation of 
Pf(t, B) to the exponential is large, meaning a large 7. For the type //regions, the unusual 
returns that happen inside the interval are associated to UPOs outside it. Such behavior is 

n 

a consequence of the existence of homoclinic tangencies [13[ inside the region. For the typ e 



III regions, the unusual returns are associated with lower-period UPOs inside them. 



13| 



When there is an UPO with low period P inside the box, p{r = P) is no longer exactly 
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equal to pnr [given by Eq. (1A3j) ] due to the fact that linearization around this low-period 
UPO does not provide the measure due to this UPO inside the box. 



III. ESTIMATION OF a, (3*, r™ P n , u, AND 7 

To estimate a and /3*, we make 7=0 (treating the power-law term as a perturbation), 
and place Pf(j, B) into Eqs. ([2D and (TjJ. We arrive that 

r - & (si 

where p e {B) = jj^Zt^ p(t, B)dr and (r e ) defined as in Eq. (j7|). 

As discussed in [13] , the exponential part of the distribution reflects the stochastic 
nature of the returns associated with UPOs inside B. It is interesting to know how 
much the exponential term P* exp [— ar'] inside pf(t,B) deviates from the Poisson function 
Y^-exp [In (1 — p)r'}, a function that describes the probability with which returns happen 
if one considers trajectories generated by uncorrelated processes. For that, we calculate fix 
and p 2 such that 

'' 1 exp [In (1 - p 2 )r'] ^ (3* exp [-ar'] (10) 



1 - fii 



We arrive that p\ = and p,i = a. The continuous part of Pf(t > t tjpo) deviates 
from the Poisson distribution (for which pi=p 2 ) whenever pi 7^ p 2 - Such deviation becomes 
larger, the larger \p — a\ (= \y- — -£r\\, where p e {B) = p{B)p e {B)). The larger \p — a\ is, 
the larger T^ P n — r min (notice that p = a when r m i n = T^ P n , what happens when e — ► 0) 
and since a depends on T^ p n it is reasonable to consider that the larger (the smaller) 



UJ 



( T UPO T min) (H) 



is, the more correlated and deterministic (uncorrelated and stochastic) the returns are. 

There are three approaches to estimate T^ p n . One is by just inspecting when p(r, B) 
presents a continuous exponential decay. The other is by fitting p(r, B) considering an 
exponential function of the type p' F (r,B) = j3*e~ aT , and Ty P n is the value for which 
[p' f (t,B) — p(r, £>)] becomes smaller than some N/L (N e N, and L is the number of 
PRT observed). The last approach is the one that will be considered in this work due to 
its experimental orientation. For that one has just to notice that as r becomes large, a as 
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estimated by Eq. approaches an asymptotic value close to the value obtained by fitting 
p(r, B) using the exponential function p' F (r, B). 

To estimate 7, we consider that T max — t^Pq >> 1, and use that [1 — (-py) 7 ] — 7 In {^-), 
a « /3*, and c = 1 (so in the limit of e -> 0, p F (r', B) -> /?*e" ar '). Then, from Eq. (J3J we 
arrive that 

7 oc - — (12) 
log((r))(r) 2 

using that max (t^J = x maa . >> 1. By making a ~ (3*, we assume that the deviation of /?p 
from the exponential law (ae~ aT ) is exclusively provided by the term (3r. In other words, 
the dynamical character of the system is provided by (3r. Therefore, manifestations of the 
deterministic behavior are more evident when -X-r is maximal and that happens for when 
t' = 1, what also means when r = r m j n . 

To simplify the presentation of our results, we rewrite Eq. fll2p as 



7 cx (r) e (13) 

IV. KOLMOGOROV-SINAI ENTROPY 

Now, we describe how the Kolmogorov- Sinai (KS) entropy can be written in terms of the 
PDF. Then, we apply this formalism to a typical physical situation when UPOs cannot be 
calculated and information about them is unknown. 

Firstly, remind that by the Pesin's equality [if]], the KS-entropy equals the sum of the 
positive Lyapunov exponents. 

For uniformly hyperbolic systems jfj (see Appendix), it is valid to write that 

where Nnr(t, B) and N(t) represent the number of non-recurrent UPOs with period r 
inside B and N(r) the total number of different UPOs of period r embedded in the chaotic 
attractor. A non-recurrent UPO inside the region B is an UPO that visits this region only 
once. 

Also, for this class of systems, we can write that 

N(t) = Cexp TXHKS (15) 
where Hks represents the Kolmogorov-Sinai entropy and C is a positive constant. 



Then, placing Eq. ffT4|) in (fT5|) . we have that Hks 
of t — ► oo. Then, for a finite r we have that 



-? + ±lof 



N nr {t,B) 
P(r,8) 



in the limit 



N NR (r,B) 



(16) 
1. Even if 



H KS (t,B)< -log 

T [ p(T,B) 

As shown in Appendix lAj it is reasonable to consider in Eq. ( |T6|) that Nnr 
there are more than one non-recurrent UPO inside the region where the returns are being 
measured, making Nnr = 1 in Eq. (|T6|) will not largely affect the estimation provided 
by this equation. So, to estimate the Lyapunov exponent of experimental systems, as the 
Chua's circuit, we will consider that 



Hks(t, B)<- log 

T 



(17) 



.p(t,B). 

This equation offers a way to estimate the KS-entropy of a chaotic system without having 
to calculate UPOs, a very difficult task. 

Equation ( |T6|) is valid for uniformly hyperbolic systems of the class for which local quan- 
tities provide good approximations for global quantities (Tent map, for example). For the 
other types of uniformly hyperbolic systems and the non-uniformly hyperbolic systems (as 
the Logistic map, Henon map, and Chua's circuit), the inequality in Eq. ([TBI) should approx- 
imate well the divergence on initial conditions for trajectories departing from the considered 
region. Therefore, to obtain a good estimation of the Lyapunov exponent one should consider 
an average of this quantity taken over many regions in phase space as in 



(H K s(T))<jjrUo g 



L — r 

k=X 



N nr (t, B) 



(18) 



p{r,B k ) 

where we are considering that this average is taken over Bk regions in phase space, with 
k = l,...,L. 

Similarly, for experimental systems, we use 

L .i r i i 

(19) 

Let us now compare our result in Eq. ffTBT) with a rigorous result [l^] valid for a chaotic 
uniformly hyperbolic map (piecewise monotone maps of an interval with a finite number of 
branches and with bounded derivative of p-bounded variation with an invariant measures 
with positive entropy) that presents one positive Lyapunov exponent. For e — > 0, we have 



(Hks(t)) <y£;;log 



L * — ' r 

k=l 
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that e ~ H>{B). From Kac's lemma we have that -rpr = fi(B) and then e ~ Also, for a 
sufficiently small interval, p(r min , B) ~ /i(£>) (assuming that p(r, £>) is well described by an 
exponential distribution), then Eq. (fl71) can be rewritten as A = ~ ln ^ , which agrees with 
the result derived in Ref. [17| that A = ~ , for r min — > oo. 



V. CORRELATION FUNCTION FOR SHORT-TERM RETURNS 

The distribution p(r, B) can be used to calculate the correlation function, a quantity that 
measures how fast a system looses information about the past as it evolves. Calling <p{x) to 
be some observable measured at the position x, the correlation between the observable <fi at 
a time T [i.e., (f)(F T (x))} with the initial observation 4>(x) is given by [l| 



C((f),T) = I <f){x)(f){F T {x))dfi (20) 
(x)dfM. / 4>(F T (x))dfj, 



where d\x stands as usually to a(x)dx and cr(x) represents the invariant density from which 
the invariant measure //(£>) can be calculated by /i(£>) = f xeB a(x)dx 

In this work, we are mainly interested in understanding the behavior of the PRTs to 
regions B. So, instead of averaging the correlation of trajectories over the whole space x, we 
calculate the correlation of trajectories in B. Employing similar ideas to the ones in Refs. 
0, Q, Q, and writing the observable to be /x(£>), we have that 

C(r, B) = y[B n F-t(B)} - h{B]h{F-t{B)] (21) 

where C(r, B) measures the correlation between trajectories that visit the region B and that 
return to it after r iterations. This function is also known as the speed of mixing. 
For an invariant measure, we have that fi[B] = [i[F~ T (B)}, and thus, 

C{r,B) = n[BnF- T (B)]- fi[B} 2 . (22) 

The quantity fi[BnF~ T (B)] represents the probability measure of having trajectory points 
leaving B and returning to it after a series of returns r^, with i — 1, . . . , I such that Yl\ T i = r 
and n < r. We can write the set [Bf]F- T (B)} as a union of two sets S'US* with S'nS* = 0, 
where S' as defined in Eq. (1A5j) and S* defined as the set of points that are mapped to B 
after r iterations but with the additional fact that for each x* G S* it exists r* < r for 
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which F T * (x*) G B. In other words, S* represents the set of points that are mapped back 
to B after r iterations, excluding all the points that firstly return to B after r iterations. 

If r < 2r min , then, /i[B n F T (B)} = fi(S'). Using Eqs. and (jA51) . we arrive to that 

p[i3 fl F r (i3)] = p(£>)p(r, B), since there cannot be trajectory points that return more than 
once within this time interval. 

Then, for r < 2r min 

C(r,B)=p(B)p(r,B)-p[B] 2 (23) 
For r > 2r m j„, we can always write that 

C(r,£)<p(£)p(r,i3)-p[£] 2 , (24) 

since fi[B n F T (23)] = p(r, £)//(£) + p[S*], Notice that as r grows, fi[S*] -> /i R in Eq. (!A2l) . 

Therefore, if the density of the returns has an exponential decay with respect to time, 
so will the correlation behave, a statement that was rigorously proved for some classes of 
systems and observables in Ref. ri|. The advantage of the approach proposed here is the 
simplicity with which we can understand such a rigorous result. 

From Eqs. (123]) and (jUj) we can clearly see that the larger (p(r, B) — fi[B]) is, the slower 
the decay of the correlation function. That is exactly the case for type II and III non- 
hyperbolic regions, for r < 2r min , when the distribution p(r, B) receives a large contribution 
of the power-law term (3r. For the type I hyperbolic regions, C(r, B) decays much faster to 0, 
since p(r, B) ~ p(B). And typically, the bigger u in Eq. (TIT!) is, the bigger (p(r, B) — p[i3]). 
That is so because the smaller r min is, the larger p{r min) B). 

In Ref. f|, a similar derivation of the correlation function was proposed, but considering 
the correlation between trajectories departing from S' and arriving to B. In other words, 
the correlation between trajectories that produce first returns and that return only once to 
B. In Eqs. (1231) and ([21]), one can estimate the correlation between all trajectories departing 
from B and arriving in B including also trajectories that return more than once. 

VI. DETERMINISM AND STOCHASTICISM IN THE HENON MAP AND THE 
EXPERIMENTAL CHUA'S CIRCUIT 



The Henon's map is given by x n+ i = a — x 2 n + by n , and y n +i = x n , wit 
b=0.3, and the considered experimental Chua's circuit can be seen in Ref. 

12 



a = 1.4 and 



19|. We focus 
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FIG. 2: [Color online] Results from the experimental Chua's circuit. (A) Probability measure with 
which X n visits intervals of length e = 0.005^. The arrows (from left to right) indicate the intervals 
B±, B3, and B%, respectively. (B) First returning map X n x X n+ \ in full gray squares, pre-images 
of trajectory points located at the maximum of the returning map in stars, and the UPOs of up 
to period 6 in pluses. 

our analyzes in three regions B\ (type I), £> 2 (type II), and B 3 (type III). 

For the Henon attractor, the regions B are boxes of equal sides 2e. The region B\ rep- 
resents a box centered at (x,y) = (1.11807,0.14719), B2 a box centered at the primary 
homoclinic tangent (x,y) = (1.780098,0.09495), and B3 a box centered at a period-2 UPO 
(x,y) = (1.36612008,-0.666120078). 

For the Chua's circuit, we reconstruct the attractor using the same techniques of Ref. []J 
and make the trajectory discrete, producing a discrete time series represented by X n , which 
represents the voltage in the capacitor C\, whenever the voltage in the capacitor C2 reaches 
zero. The attractor is of the Rossler-type. In Fig. 12(A) we show the invariant measure of X n 
and in [2(B), the first returning map X n x X n+1 , by the gray full squares. The (blue) stars 
represent trajectory points located at the maximum of the return map and its pre-images. 
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FIG. 3: [Color online] (A) and (B) are results for the Henon attractor and (C) and (D) results for 
the experimental Chua's circuit. Estimation of the parameter a and r™pQ for the Henon attractor 

(A) and for the Chua's circuit (C). In (B) and (D), we show p(r,B), the fitting of p(r,B) by a 
function of the type Pf(t, B) [Eq. ([5])], and the fitting of p(r, B) by an exponential function of the 
type p' f (t,B) [Eq. ([T])] . The fittings are made by considering the transformed time r' [Eq. ([6])] 
but, in these figures, we re-transform the time back to r in order to plot the fittings together with 
the distribution p(r,B). 

Every pre-image is located at a X n point for which the probability density is large in (A), 
what typically happens for non-uniformly hyperbolic systems. The plus symbol indicates 
the place of the lower-period UPOs (up to period 6). In this circuit, B represents intervals 
of size e. B\ represents the interval centered at X n = —0.64, 23 2 the interval centered at 
X n = 0.05, and 3% the interval centered at the position of a period-2 UPO X n = —0.586637. 
These points where the intervals are positioned are indicated by the arrows in Fig. [2(A). 

We now make a detailed analysis of the PRTs to the box £>2, for the Henon map [in Figs. 
EKA-B)], and the PRTs to the interval B 2 , for the Chua's circuit [in Figs. E](C-D)]. In (A) 
and (C), we estimate the value of a as we consider different values for T™ p n , in Eq. We 
consider the smallest value of t™™ for which the value of a " converges" . The arrow indicates 
these values. For the Henon attractor, t™ P q = 28 and for the Chua's circuit, T™ % p Q = 13. In 

(B) and (D), we show p(r, B), Pf(t, B), and a fitting of p(r, B) by an exponential function of 
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FIG. 4: The fittings of 7 (vertical axis) versus (r) in a log-log graph, for the distributions pp 
of Henon map in (A-C) (left column) and for the experimental Chua's circuit in (D-F) (right 
column), in the form 7 oc (r) e [see Eq. (I13p ], The value of the power 6 is shown inside the figures. 
Experimentally, the value of (r) is over-estimated, which leads to a smaller 9 exponent. 

the type of p' F (r, B) [Eq. (pQ)]. For the Henon attractor (B), using the distribution p F (r, B) in 
the integral of Eq. Q produces the value 0.994826 while using p' F (r, B) produces 0.987701. 
For the Chua's circuit (D), using the distribution p F (r, B) in the integral of Eq. ([3]) produces 
the value 0.78 while using p' F (r, B) produces 0.68. Reminding that the integrals in Eq. (J3j) 
should provide 1, it is clear that the proposed distribution in Eq. ([5]) fits better p(r, B). 

In Fig. HI we show the relation between 7 and (r) for S 1; B2, and B3, in the Henon 
attractor [first column, (A-C)] and in the Chua's circuit [second column, (D-F)]. We have 
numerically obtained that for the Henon map, 7 oc (r) e with 9 = —0.6 [excluding the region 
£>3, whose results are shown in (C)], and for the Chua's circuit 7 oc (t) 9 with 9 < —0.9. 
Remind that the larger is \9\, the smaller /3r. 

For type II non-hyperbolic regions, as e is decreased, typically we expect that the unstable 
manifold of the UPO outside the region will still belong to the region, which leads to a small 
\9\ value. 

Table fl] shows estimates for the positive Lyapunov exponent using the right-hand side 
of Eq. ( fl~6l) for the Henon map and Eq. (JTTJ) for the Chua's circuit. B\ represents a 
hyperbolic type I region and B2 a non-hyperbolic type III region, with different sizes e. 
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TABLE I: Estimation of the Lyapunov exponent. For the Henon map, we typically find that the 
UPO with the smallest period has a period larger than the observed r m j n , a consequence of the 
non-hyperbolic character of this map [see Fig. [6^B) and[6^D)]. Thus, for calculating the Lyapunov 
exponent, we consider in Eq. (fl~6j) r = P m in, where P m in is the lowest period of all UPOs inside 
B. For all regions studied in this map, the number of non-recurrent UPOs is within the interval 
Nnr = [1; 2] as expected. For the experimental Chua's circuit, we calculate the Lyapunov exponent 
considering in Eq. (|17p r equal to the first Poincare return for which the PDF presents its third 



maximum. We also assume that N^r = 1. The 
is also calculated using the technique in Ref. 



positive Lyapunov exponent of the Henon attractor 



the largest Lyapunov exponent of the Chua's 



20 ] , with the code of the Tisean package 



2j. 



161 ] anc 

circuit is also calculated using the technique in Ref. 
We obtain that the positive Lyapunov exponent of the Henon map is 0.419/iteration and the one 
for the Chua's circuit is 0.52/cycle. The value of r used is the number between parentheses. 



e - Henon 


0.02 


0.01 


0.005 


0.001 


0.0008 


Bi 


0.32(21) 


0.40(23) 


0.40(23) 


0.39(23) 


0.39(23) 


B 2 


0.61(18) 


0.33(18) 


0.36(18) 


0.40(22) 


0.29(29) 


e - Chua 


0.1700 


0.1450 


0.1200 


0.0950 


0.0700 


Bi 


0.43(7) 


0.42(7) 


0.42(7) 


0.42(7) 


0.46(9) 


B 2 


0.58(5) 


0.39(7) 


0.58(5) 


0.25(13) 


0.27(13) 



For the Henon attractor, a numerical calculation of this exponent provides A=0.419 16], 
and for the experimental Chua's circuit, A=0.52/cycle 20, 12JJ- The unit of the Lyapunov 
exponent for this circuit is in [1/cycles]. That is due to the fact that for the calculation 
of this exponent we have used the discrete time series X n , which represents the voltage in 
capacitor C\ whenever the voltage of capacitor C 2 is zero. 

In both systems, the Lyapunov exponents estimated from Eqs. (fl~6l) [for the Henon map] 
and (fTTl) [for the Chua's circuit] using returns measured in the non- hyperbolic regions B 2 
produce worse estimates than the ones measured in the hyperbolic regions B%. That is to be 
expected since our estimations are valid for uniformly hyperbolic systems. For both systems 
and all regions, the estimation of the Lyapunov exponents produce better results as the size 
of the regions decrease and T min becomes large. That is also to be expected since as the size 
of the box decreases the chance that a region has a hyperbolic character increases. 
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Since both the Henon map and the Chua's circuit are non-uniformly hyperbolic systems, 
in order to obtain better estimates for the Lyapunov exponents we need to calculate an 
average value considering many regions in phase space. 

To firstly illustrate how an average value provides a better estimation for the real value 
of the Lyapunov exponent, we average all the values shown in Table [B For the Henon map, 
we obtain that A = 0.39 (the real Lyapunov exponent is 0.419) and for the Chua's circuit 
we obtain that A = 0.42 (the real value is 0.52). 

Then, we calculate (H KS ) for the Henon map, using Eq. (fl8l) and considering L =500 
regions with e=0.005 and r = P m i n , where P m in is the lowest period of all UPOs inside each 
region. We obtain that (H^s) — 0.464 (compare with the real value 0.419). These L =500 
regions are centered in points of a 500 long trajectory. Using Eq. (fT9~|) for the same previous 
conditions, we obtain (Hks) = 0.441. For the Chua's circuit, we calculate (Hks) using 
Eq. (TT9T) considering L =25 regions B with e=0.067 and r equal to the second largest first 
Poincare return. We obtain that (Hks) — 0.42 ±0.1 (compare with the real value 0.52). 

In Fig. El we show the decay of the correlation function with respect to r. Notice that 
this function is local and reflects the decay of correlation of PRTs to a particular interval. 
As expected, the correlation function decays faster in the type I hyperbolic regions [(A) and 
(C)] than in the type II non-hyperbolic regions. In addition, in general, the larger e is, the 
larger the value of the correlation function, a consequence of the fact that the larger e is, 
the larger both 9 and u are. 

VII. CONCLUSIONS 

We have studied the Poincare recurrence time (PRT) of chaotic systems, a quantity that 
measures the time interval between two consecutive visits of a trajectory to a finite size 
region in phase space. The motivation in studying PRT is that for some systems the only 
possible measure one can make is the time interval between two events. But one still wants 
to understand what kind of dynamics is behind the generation of these returns. 

If the region in phase space has an arbitrarily small size, for systems that have strong 
mixing characteristics, a long return time contains no longer information about the future 
returns, which means that the consecutive series of returns lose their correlation and they 
behave as if they had been generated by a completely random process. At such a situation, 
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FIG. 5: Correlation function with respect to r for the Henon map (A-B) and the experimental 
Chua's circuit (C-D). In (A) and (C), we consider the type I hyperbolic regions Bi, and in (B) and 
(D), we consider the type II non- hyperbolic regions 82- 



the probability distribution of the PRTs approaches a Poisson distribution [15|, |22| and few 
can be said about the dynamical manifestations of the data by only considering the form of 
this distribution. 

Our approach is mainly devoted to characterize a chaotic system (obtaining relevant 
invariant quantities) considering regions with finite size and short return times. We show 
how to calculate the short-term correlation function [see Eq. (|23l) ] and the Kolmogorov-Sinai 
entropy [see Eqs. f|T7|) and ffl9l) ]. using the distribution of PRTs, a quantity that can be 
easily accessible in experiments. 

As the region in phase space where the returns are being measured becomes larger, the 
first Poincare returns reflect more the deterministic nature of the system. That leads to a 
large deviation of the density of the returns from the Poisson [see Eq. 03)] an d a slower 
decay of correlations for short-term returns [see Eq. fl23|) and Fig. [5]. 

These characteristics can be advantageously used to characterize the level of determinism 
in chaotic systems. More specifically, the larger 9 is in Eq. ffTTl) . the more deterministic 
is the considered chaotic system under the point of view of the PRTs. Another quantity 
is to = T™ l p — r min , where T^ P n is the return value for which the probability distribution 
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of the PRTs presents a continuous decay with the increasing of the return time, and r m j n 
is the return with the minimal length. The larger u is, the larger is the time span within 
which discontinuities are observed in the distribution of returns. Each discontinuity identifies 
particular returns, which makes them to be highly predictable. As shown in Sec. |V] these 
previous intuitive ideas are indeed correct, i.e., larger 9 and uj result in a slower decay of 
correlations. 

The strategy of using a value of r that is larger than r m j„ in Eqs. ( Tl8|) and (1T91) is due 
to the fact that in non-hyperbolic systems as the Henon map and the experimental Chua's 
circuit the density of the PRTs for r = r min , i.e. p(T m i n ,B), receives a contribution coming 
from UPOs outside of B, which violates the conditions under which these equations are 
derived. When the UPOs are known as it is the case for the Henon map, we use r = P m i n , 
in Eqs. ( ITS]) and (JT9"j) . where P m in represents the period of the UPO inside B with the lowest 
period. When the UPOs are unknown as it is the case for the experimental Chua's circuit, 
we use t > T min . 

Our approach can be straightforward applied to the treatment of data coming from 
complex systems, and the reason lies on the formulas that we have used. As one can check 
the 3 most relevant formulas used in our work [Eqs. (]lip . (|19p . and (1231) ] depend only on 
the probability of the Poincare return times to finite size intervals with short time-length, a 
quantity that can be easily and quickly accessible from measurements and that, in principle, 
does not require the existence of chaos. 



For some technological applications as in the control of chaos (23j, one does not need 
to know with high precision the position of an UPO, but rather its unstable eigenvalue. 
Equation (1A8I) offers a trivial way to calculate this quantity by measuring the probability 
with which the fastest first Poincare return to a sufficiently small interval happens. Naturally, 
this equations also offers a simple way to obtain estimations for the local first derivative of 
a system, a quantity often needed to make local models of complex systems. If the system 
is higher- dimensional then L m , n in Eq. (1A8I) should refer to the product of all the unstab 



eigenvalues of a single non-recurrent UPO appearing in the region, as explained in Ref. 24]. 
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APPENDIX A: THE DENSITY OF RETURNS p AND THE NON-RECURRENT 
UPOS 



In Ref. 24| it is derived a formula that relates the probability measure of a D-dimensional 
box B with the unstable eigenvalues of the UPOs inside it. More specifically, 



KB) 



lim 

P^oo 



L J 



(Al) 



where Lj(P) = Lij(P)L 2 j(P) . . . Ljjj(P), and L u j(P) (u = 1, . . . , U) represent the U posi- 
tive eigenvalues (larger than 1) of the j fixed point located in B of the P-fold iterate of the 
map represented by F p (i.e., the fixed points are period-P UPOs that belongs to B), where 
P tends to infinity. In a general situation, there are U positive eigenvalues and the box 
is D-dimensional. In the following, we consider D=2, and U=l (there is only one positive 
Lyapunov exponent). 

Equation (lAlj) was demonstrated to hold for mixing hyperbolic (axiom A) attractors and 
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25|], for 



was shown numerically to hold for the non-hyperbolic Henon attractor in Refs. 
UPOs of moderately large period P = 30. 

As done in Ref. [6j , we rewrite the right-hand side of Eq. (jAlft as a sum of two terms, 
without taking the limit of P — > oo but for a finite P: 



(A2) 



where y, R = J2j[Lf(P)] 1 and 



(A3) 



where Lf(P) are the unstable eigenvalues of the so called recurrent UPOs that return to 
B more than once before completing their cycles, and L^ R {P) are the unstable eigenvalues 
of the so called non-recurrent UPOs that return to B only once. So, while fi^R measures 
the contribution to the measure due to chaotic trajectories associated with non-recurrent 
UPOs, hr measures the contribution to the measure due to chaotic trajectories associated 
with recurrent UPOs. 
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As shown in Ref. Q], there is a clever way to relate the density p(r, B) with the measure 
of the attractor associated with UPOs inside B by 

p(r,B)=fi NR (r,B) (A4) 

The term Pnr(t, B) can be represented in terms of space averages by 

^ NR{T ' B) = W) (A5) 

where fi(S') represents the measure of the set S' (the part of the measure inside B due to 
the set 5") and 5" represents the set of points of the attractor that returns firstly to B after 
r iterations. More rigorously, representing by F the transformation that generates a chaotic 
set A and given a subset B C A, then S' G A and S' = F~ T {B) n B such that there is not 
t* <r for which F T * (B)nB^ 0. 

But the right-hand side of Eq. flA5j) can be estimated by Nn ^^ , and therefore, 

/ m N nr (t,B) 
Vnr(t, B) = v , (A6) 

where N^ r (t, B) and N(t) represent the number of non-recurrent UPOs with period r 
inside B and N(r) the total number of different UPOs of period r embedded in the chaotic 
attractor. 

Now, using Eq. (1A3j) . it is clear that 

P(r, B) > (A7) 

where L min represents the unstable eigenvalue with the lowest amplitude within all UPOs 
with period r. 

In Fig. [6J we show the values of L m j n (r) _1 and the density of returns p in the Henon 
attractor, for the regions B\ [in (A)], i5 2 [in (B)], B3 [in (C)], and also for a region centered 
at a period-8 UPO [in (D)] at the position (x, y)=(1.496703, -0.545333), denoted by B 4 . 

We see that for r ~ T min , inequality in Eq. (1A7|) is close to an equality and we can write 
that 

The reason is that for such a case, there is only a few (a number of the order of 1) 
non-recurrent UPOs with period r m i n . It is easy to understand that by using the following 
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FIG. 6: By empty circles we show 1/L m j n (r) and (red) crosses the value of p(r, £>), for the Henon 
attractor for different regions. B\ [in (A)], B2 [in (B)], and B3 [in (C)], and also for a region centered 
at a period-8 UPO [in (D)]. e = 0.02 in all figures. The arrow points to the value of r m j n . 

argument. Imagine a sufficiently small region centered around an UPO with period P = 1. 
Clearly r m j„ = 1 and there will be only one non-recurrent UPO with period P—l. Now, 
consider a region around an UPO with period P=2. Similarly, r min = 2 and there will be 
only one UPO with period 2. Typically, for sufficiently small regions, there will be only one 
non- recurrent UPO inside the regions if P ~ T min . 

Provided that the UPO is hyperbolic, the uniqueness of the UPO in the small region is 



garanteed by the Hartman-Grobman Theorem [26j and the size of the region in which the 
uniqueness of the UPO can be garanteed is related to the strength of the hyperbolicity of 
the UPO. 

To illustrate that, we consider the Logistic map [x n+ i = bx n (l — x n )], whose bifurcation 
diagram is shown in Fig. [7(A), constructed considering 100 parameter values b within 
the parameter range [3.6,3.99]. In Fig. UTB), we show the number of non-recurrent UPOs, 
denoted by N NR , for UPOs with period P = r min and intervals with size e = 0.001, randomly 
selected such that r min 6 [10,14]. For the large majority of intervals, N NR = 1. Finally, 
in [7(C), we show the values of (Hks) calculated from Eq. (fT9l) for this parameter range. 
Notice that despite the choice of N NR = 1, the value of (H KS ) is a good estimation for the 
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positive Lyapunov exponent of this map, indicated in this figure by A. 
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FIG. 7: [Color online] Results for the Logistic map: x n+ \ = bx n (l — x n ). (A) Bifurcation diagram 
of the Logistic map (green points). Empty squares represent intervals of size e = 0.001 randomly 
selected such that r m j n G [10,14]. We consider 100 parameter values within the range [3.6,3.99] 
and for each parameter we consider one interval. (B) the number Nnr of non-recurrent UPOs with 
period P = T m i n inside each one of the 100 intervals. (C) The value of the Lyapunov exponent 
(black tick line), denoted by A, and (Hks) (red pluses) calculated from Eq. ([19]) . 
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